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2 Abstract 
O 

£N) A 3-dimensional non-commutative oscillator with no mass term but with a certain momentum- 

dependent potential admits a conserved Runge-Lenz vector, derived from the dual description in 
momentum space. The latter corresponds to a Dirac monopole with a fine-tuned inverse-square 
plus Newtonian potential, introduced by Mcintosh, Cisneros, and by Zwanziger some time ago. 
The trajectories are (arcs of) ellipses, which, in the commutative limit, reduce to the circular 

-t— > 

hodographs of the Kepler problem. The dynamical symmetry allows for an algebraic determination 

<D 

°f the bound-state spectrum and actually extends to the conformal algebra o(4, 2). 
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The "Berry phase" induced in the semiclassical analysis of a Bloch electron in a 3- 
dimensional crystal makes the coordinates non-commutative pQ, 

{xi,Xj} = e ijk Qk, {xi,Pj} = {Pi,Pj} = 0, (1) 

where 0j = Qi(p). Then the Jacobi identity requires the consistency condition V^- = 1 . 
Choosing, for example, the non-commutative vector aligned in the third direction, 0, = 
98is, 9 = const ., the 3D theory reduces to the planar mechanics based on "exotic" Galilean 
symmetry [2]. For the Kepler potential V oc r~ 1 , for example, interesting results, including 
perihelion precession, can be derived [3]. Other applications concern the Quantum Hall 
Effect 

Such a choice only allows for axial symmetry, though. Full rotational symmetry can, 
however, be restored by choosing instead to be a "monopole in ^-space" [5] , 

Qi = 6 % , 9 = const . , (2) 
pi 

where p = \p\ . We record for further use that the associated symplectic structure has an 
extra, "monopole" term, Q = dpi A dxi + (9/2p 3 )eijkPidpj A dp k , that of a mass-zero, spin-6 1 
orbit of the Poincare [and indeed of the o(4, 2) conformal] group [6]. 

In this Letter we study the 3D mechanics with non-commutativity ([2]), augmented with 
a Hamiltonian of the form H = p 2 /2 + V(x, p). Note that the potential may also depend on 
the momentum, p. The equations of motion read 

M dV Pk dV . dV 

Xi = pi + — + 9e ijk — — , Pi = -—. (3) 

OPi p 6 OXj OXi 

Note, in the first relation, the "anomalous velocity terms" due to our assumptions. 

We are particularly interested in finding conserved quantities. This task is conveniently 
achieved by using van Holten's covariant framework [T]. which amounts to searching for 
an expansion in integer powers of the momentum, Q = Cq(x) + Ci(x)pi + ^Cij(x)piPj + 
^Cijk(x)piPjPk+ ■ ■ ■ Requiring Q to Poisson-commute with the Hamiltonian yields an infinite 
series of constraints. The expansion can be truncated at a finite order n, however, provided 
the Killing equation is satisfied, %C?i J ,.,i B ) = 0, when we can set Cj 1 ...i n+1 ... = 0. 

1 For a more general theory which also includes magnetic fields, see, e.g., pfl[2]- For simplicity, the mass 
has been chosen unity. 
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Let us assume that the potential has the form V(r,p), and try to find the conserved 
angular momentum, associated with the Killing vector C = n x x, which represents a space 
rotation around n. An easy calculation shows, however, that the procedure fails to work, 
owing to the "monopole" term, which is not an integer power of p. We propose therefore to 
work instead in "dual" (i.e. momentum) space, and search for conserved quantities expanded 
rather into powers of the position, 

Q = C (p) + Ci(p)Xi + yCij(p)XiXj + yCijkipjXiXjXk ... (4) 

Then van Holten's algorithm is replaced by 

Ci[Pi + ^) = o(0) 



dpi/ 

de ijk^Ci - — + dj A Pi + ^— =0 0(1 

r Or \ p 6 opj J \ dpi J 

— 5~ e ^\ e ijmC ik + eikmCij) - (-5— + -^-) ) +C ijk [Pi + — ) =0 o(2 

r Or \ p A v ' v opj Op k ' J \ dpi J 

IdV f n 

7\ V^itlim^ljk + Cljm^lki + ^lim^lij) ~ { ~^ 1 r ~~ ^ ) + 

r or V p opk opi opj J 

dV\ 



Cujk [Pi+Qp-J = o(3) 



where r = |x|. Then, for the dual Killing vector C = n x p, the algorithm provides us with 
the conserved angular momentum, 

J = L — 9p = xxp — Op P = p/p, (5) 

which is what one would expect, due to the "monopole in p-space", whereas the NC param- 
eter, 9, behaves as the "monopole charge" [8]. 

The next step is to inquire about second order conserved quantities. The usual Runge- 
Lenz vector is generated by the Killing tensor = 25^ where n is some 

fixed unit vector [7]. Not surprisingly, the original procedure fails once again. The dual 
procedure Q works, though. The two-tensor 

Cij = 26ijn ■ p - riipj - rijPi, (6) 

verifies the dual Killing equation of order 3. Then the order-2 equation yields Ck = 9 ^ nx ^ >k . 
Inserting into the first-order constraint and assuming d r V 7^ 0, the constraint is satisfied 



with C = an ■ p, provided the potential and the Hamiltonian are 



V 




H 



T + 2pi 



+ 



a 




(7) 



respectively, where a is an arbitrary constant. Then the dual algorithm provides us with 
the Runge-Lenz-type vector 

K = x x J — ap. (8) 



Its conservation can also be checked by a direct calculation, using the equations of motion, 



Note that the —p 2 /2 term in the potential cancels the usual kinetic term, and our system 
describes a non-relativistic, non- commutative particle with no mass term in an oscillator 
field, plus some momentum- dependent interaction. Writing the Hamiltonian as H = x 2 /2 + 



when p = —8 2 /a, which plainly requires a < 0. 

It is easy to understand the reason why our modified algorithm did work : calling p 
"position" and — x "momentum", the system can also be interpreted as an "ordinary" (i.e. 
massive and commutative) non-relativistic charged particle in the field of a Dirac monopole 
of strength 9, augmented with an inverse-square plus a Newtonian potential. This is the 
well-known "Mcintosh- Cisneros - Zwanziger" (MICZ) system [6j [9], for which the fine- 
tuned inverse-square potential is known to cancel the effect of the monopole, allowing for a 
Kepler- type dynamical symmetry [HI E] ■ The angular momentum, (|5]) , and the Runge-Lenz 
vector, (18]), are, in particular, that of the MICZ problem [9] in "dual" [momentum] space. 

The conserved quantities provide us with valuable information on the motion. Mimicking 
what is done in the MICZ case, we note that J ■ p = —6 implies that the vector p moves 
on a cone of opening angle arccos(— 9/ J). On the other hand, for iV = a J — 9K we have 
N ■ p = 9(J 2 — 6 2 ) = 9L 2 a constant, so that the p-motion is in the plane perpendicular to 
N. The trajectory in p-space belongs therefore to a conic section. 

For the MICZ problem, this is the main result, but for us here our main interest lies in 
finding the real space trajectories, x(t). By ([9]), this amounts to finding the [momentum-] 
"hodograph" of the MICZ problem. Curiously, while the hodograph of the Kepler problem 
is well-known [it is a circle or a circular arc], we could not find the corresponding result in 




(9) 



(6 2 /2)(p- 1 +a/e 2 ) 2 



a 2 /28 2 shows, moreover, that H > —a 2 /26 2 with equality only attained 
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the (vast) literature of MICZ. Returning to our notations, we note that due to N ■ x — 0, 
x(t) also belongs to the same oblique plane, whose normal is N. 

The problem is conveniently studied in an adapted coordinate system. One proves indeed 
that 

|?J,jfe} = {j^K x J, j^-(29HJ+aK), j^{af-9K)} (10) 

[where A 2 = a 2 + 2H9 2 , e 2 = a 2 + 2HJ 2 and L 2 = J 2 — 9 2 ] is an orthonormal basis. (In k 
we recognize, in particular, N/N.) 

Firstly, projecting onto these axis, p z = p '■ k = 9L/\X\ = const, while for p x = p '■ i and 
Py = P ■ J we find 

(Py + w) vl = fn 
\2/4 H 2 L 2/2H ' 1 ; 

which is the equation of a hyperbola or of an ellipse, depending on the sign of H . (For 

vanishing H one gets a parabola). This confirms what is known for the MICZ problem [9]. 

Next, for 

X = £-i=--t-±(H-—), Y = £.j=-Ll—?. (12) 
|e| 2p \e\ p 



an easy calculation yields 



lelL.o a 2 L 2 n L 2 a 2 . „. 



[completed with Z = x ■ k = 0] which is an ellipse, since A 2 = a 2 + 2H6 2 > 0. The center 
has been shifted to — \e\L/J 2 along the axis i. The major axis is directed along j. 

Note that, unlike as in p-space, the x-trajectories are always bounded. When the energy 
is negative, H < [which is only possible when the Newtonian potential is attractive, a < 0], 
the whole ellipse is described. When H > 0, [which is the only possibility in the repulsive 
case a > 0] it is only the arc between the tangents drawn from the origin which is obtained. 
When the non-commutativity is turned off, 9 — > 0, the known circular hodographs of the 
Kepler problem are recovered. As a — > 0, the trajectory becomes unbounded, and follows 
the y-axis. 

So far, we only discussed classical mechanics. Quantization is now straightforward using 
the known group theoretical properties of the MICZ problem in dual space. The non- 
commutativity [alias monopole charge] 9 has to be an integer or half integer - this is indeed 
the first indication about the quantization of the NC parameter - ; wave functions should be 
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FIG. 1: When the Newtonian potential is attractive, a < 0, one can have H < 0, and the trajectories 
are full ellipses. The origin is inside the ellipse. When H > 0, which can arise either in the high- 
energy attractive or in the repulsive (a > 0) case, the origin is outside our ellipse. The particle is 
confined to an arc [denoted with the heavy line] of the ellipse, determined by the tangents drawn 
from the origin. For H = 0, the origin lies on the ellipse, and "motion" reduces to this single point. 

chosen in the momentum representation, ip\P)- The angular momentum, J, and the rescaled 
Runge-Lenz vector, K/^/2\H\, close into o(3, l)/o(4) depending on the sign of the energy. 
In the last case, the representation theory provides us with the discrete energy spectrum [in 
units h = 1] 



where n = 0, 1, ...,/ = 0, 1, ... , with degeneracy n 2 — 6 2 . The same result can plainly be 
derived directly by solving the Schrodinger equation in p-space 0. 

Moreover, the symmetry extends to the conformal o(4, 2) symmetry, due to the fact that 
the massless Poincare orbits with helicity 9 are in fact orbits of the conformal group, cf. [6]. 

In conclusion, we observe that in most approaches one studies the properties (like tra- 
jectories, symmetries, etc.) of some given physical system. Here we followed in the reverse 
direction: after positing the fundamental commutation relations, we were looking for poten- 
tials with nice properties. This leads us to the momentum-dependent potentials ([7]), which 
is indeed a Mclntosh-Cisneros — Zwanziger system [H] in dual space. Unlike as in a constant 
electric field pi)], the motions lie in an (oblique) plane. The particle is confined to bounded 
trajectories, namely to (arcs of) ellipses. 

The best way to figure our motions is to think of them as analogs of the circular 



'n 





(14) 



hodographs of the Kepler problem to which they indeed reduce when the non-commutativity 
is turned off. For H < 0, for example, the dual motions are bound, and the velocity turns 
around the whole ellipse; for H > instead, our motion along a finite arc, starting from 
one extreme point and tending to the other one at the end of the arc, corresponds to the 
variation of the velocity in the course of a hyperbolic motion (of a comet, or in Rutherford 
scattering) in dual space. 

There is also some analogy with particles with torsion [11] . 

Our system, with monopole-type non-commutativity (|2j), has some remarkable properties. 

Momentum-dependent potentials are rather unusual in high-energy physics. They are, 
however, widely used in nuclear physics, namely in the study of heavy ion collisions, where 
they correspond to non-local interactions [12]. Remarkably, in NC field theory, a 1/p 2 
contribution to the propagator emerges from UV-IR mixing. See the middle equation in 
equation (22) page 11 of Ref. [13]. 

The absence of a mass term should not be thought of as the system being massless; it is 
rather reminiscent of "Chern-Simons dynamics" [I]. 

One can be puzzled how the system would look like in configuration space. Trying to 
eliminate the momentum from the phase-space equations (|9| in the usual way, which amounts 
to deriving x w.r.t. time and using the equations for p, fails, however, owing to the presence 
of underived p in the resulting equation. This reflects the non-local character of the system. 

One can, instead, eliminate x using the same procedure, but in dual space. This yields 
in fact the equations of motion of MICZ in momentum space, 

.. J 2 a a _ 9 -> 

P=^ + ^, P= -sP~ ~^ J - ( 15 ) 



Are these equations related to a theory with higher-order derivatives of the type [13] ? 
The answer is yes and no. The clue is that time is not a "good" parameter for Kepler-type 
problems, owing to the impossibility of expressing it from the Kepler equation [15] . This is 
also the reason for which we describe the shape of the trajectories, but we do not integrate 
the equations of motion. A "better" parameter can be found along the lines indicated 



by Souriau [TBI [17J and then, deriving w.r.t. the new parameter, transforms (15) into a 
fourth-order linear matrix differential equation, which can be solved. 

It is, however, not clear at all if these equations derive from some higher-order Lagrangian, 
and if they happen to do, what would be the physical meaning of the latter. 
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The fourth-order equations do certainly not come from one of the type in Ref. [H] : the 
latter lives in fact in two space dimensions and has constant scalar non-commutativity 9 - 
while our system is 3-dimensional and has a momentum- dependent vector Q(p), given in 

It is tempting to ask if the relation to the "closest physical theory" with a momentum 
- dependent potential, namely nuclear physics, can be further developed and if similar (su- 
per) symmetries can be found also in nuclear physics. Once again, the answer seems to be 
negative, though : while dynamical symmetries do play a role in nuclear physics [IS], those 
used so far do not seem to be of a momentum-dependent Keplerian types. 

Acknowledgments 

PAH is indebted to the Institute of Modern Physics of the Lanzhou branch of the Chinese 
Academy of Sciences for a special professorship, and PMZ is indebted to the Laboratoire de 
Mathematiques et de Physique Theorique of Tours University for hospitality. JPN would like 
to thank the Region Centre for a doctoral scholarship. We are grateful to M. Plyushchay 
and to D. O'Connor for useful correspondence. 



[1] M. C. Chang and Q. Niu, "Berry Phase, Hyperorbits, and the Hofstadter Spectrum." Phys. 

Rev. Lett. 75, 1348 (1995). For a review, see D. Xiao, M. C. Chang and Q. Niu, "Berry Phase 

Effects on Electronic Properties," |arXiv:0907!202T1 
[2] It is impossible to give here an even approximately complete list of relevant contributions. The 

original references are C. Duval and P. A. Horvathy, "The "Peierls substitution" and the exotic 



Galilei group," Phys. Lett. B 479, 284-290 (2000). |hep-th/0002233|; "Exotic galilean sym- 
metry in the non-commutative plane, and the Hall effect," Journ. Phys. A34, 10097 (2001). 
hep-th/0106089|; V. P. Nair and A. P. Polychronakos, "Quantum mechanics on the non- 



commutative plane and sphere," Phys. Lett. B 505 (2001) 267 [arXiv:hep -th/0011172|. For a 
recent review, see, e.g. P. A. Horvathy, L. Martina and P. C. Stichel, "Exotic galilean symme- 



try and non-commutative mechanics," SIGMA 6 (2010), 060 arXiv: 1002.4772 and references 
therein. 



8 



[3] J. M. Romero and J. D. Vergara, "The Kepler problem and non commutativity," Mod. Phys. 



Lett. A 18 (2003) 1673 |arXiv:hep-th/0303064i. 
[4] G. V. Dunne, R. Jackiw and C. A. Trugenberger, "Topological (Chern-Simons) Quantum 
Mechanics," Phys. Rev. D 41 (1990) 661; In the non-commutative case, see P. A. Horvathy, 



"The non-commutative Landau problem." Ann. Phys. (N. Y.) 299, 128-140 (2002) |hep- 
th/0201007 . 

[5] A. Berard and H. Mohrbach, "Monopole in momentum space in noncommutative quantum 



mechanics," Phys. Rev. D 69 (2004) 127701 arXiv:hep-th/0310167|. The analogy with anyon 



has been discussed in M. S. Plyushchay, "Monopole Chern-Simons term: Charge monopole 



system as a particle with spin." Nucl. Phys. B589 413 (2000), [hep-th /0004032] 
[6] P. A. Horvathy, B. Cordani and L. Feher : "Kepler-type dynamical symmetries of long- 
range monopole interactions." Journ. Math. Phys. 31, 202 (1990). See also M. S. Plyushchay, 
"Massless point particle with rigidity," Mod. Phys. Lett. A 4, 837 (1989). 
[7] J. W. van Holten, "Covariant hamiltonian dynamics," Phys. Rev. D 75, 025027 (2007) 
arX iv:hep-th/0612216| . Using Killing tensors has been advocated by B. Carter, "Global 
structure of the Kerr family of gravitational fields," Phys. Rev. 174, 1559 (1968). For re- 
cent applications, see e.g., P. A. Horvathy and J. P. Ngome, "Conserved quantities in a 



non-abelian monopole field," Phys. Rev. D 79 (2009) 127701 [arXiv: 0902.0273|. J. P. Ngome, 
"Curved manifolds with conserved Runge-Lenz vectors," J. Math. Phys. 50 (2009) 122901 
arXiv:0908.1204|. M. Visinescu, "Higher order first integrals of motion in a gauge covariant 



Hamiltonian framework," Mod. Phys. Lett. A 25 (2010) 341 arXiv:0910.3474|. J. P. Ngome, 



P. A. Horvathy and J. W. van Holten, "Dynamical supersymmetry of spin particle-magnetic 



field interaction," J. Phys. A: Math. Theor. A43 (2010) 285401, arXiv: 1003.0137 T. Igata, 



T. Koike and H. Ishihara, "Constants of Motion for Constrained Hamiltonian Systems: A 



Particle around a Charged Rotating Black Hole," arXiv:1005.1815 [gr-qc]. 
[8] Jose L. Cortes, Mikhail S. Plyushchay, "Anyons as spinning particles," Int. J. Mod. Phys. All 



3331 (1996). |hep-th/9505117| 
[9] H. V. Mcintosh and A. Cisneros, "Degeneracy in the presence of a magnetic monopole." Journ. 
Math. Phys. (N.Y.) 11, 896 (1970); D. Zwanziger, "Quantum field theory of particles with 
both electric and magnetic charges." Phys. Rev. 176, 1489 (1968). 
[10] P. A. Horvathy, "The anomalous Hall effect in non-commutative mechanics," Phys. Lett. A 



9 



359, 705 (2006). [cond-mat/0606472] . 

[11] M. Plyushchay, "Relativistic particle with torsion and charged particle in a constant electro- 
magnetic field : identity of evolution," Mod. Phys. Lett. A10, 1463 (1995) [hep-th/0930914| 7]. 

[12] C. Gale, G. Bertsch, and S. Das Gupta, "Heavy-ion collision theory with momentum- 
dependent interactions," Phys. Rev. C 35, 1666 (1987); C.B. Das, S.Das Gupta, Charles 
Gale, and Bao-An Li, "Momentum dependence of symmetry potential in asymmetric nu- 
clear matter for transport model calculations," Phys. Rev. C 67 (2003) 034611; "Effects of 
momentum-dependent symmetry potential on heavy-ion collisions induced by neutron-rich 
nuclei," Bao-An Li, Ch. B. Das, S. D. Gupta, Ch. Gale, Nuclear Physics A 735 (2004) 563. 

[13] S. S. Gubser and S. L. Sondhi, "Phase structure of non-commutative scalar field theories," 
Nucl. Phys. B605 395 (2001) jhep-th/0006119] . 

[14] Lukierski J., Stichel P. C. and Zakrzewski W. J., "Galilean-invariant (2 + l)-dimensional 
models with a Chern-Simons-like term and d = 2 noncommutative geometry," Annals of 
Physics 260 (1997), 224, |hep-th/9612017|; "Noncommutative planar particle dynamics with 
gauge interactions", Ann. Phys. (N.Y.) 306 (2003), 78 [hep-th/0207149] . 

[15] B. Cordani, The Kepler Problem. Group Aspects, Regularization and Quantization. With an 
Application to the Study of Perturbations. Birkhauser: Basel (2003). 

[16] Jean-Marie Souriau, "Global geometry of the two-body problem," (In French), CPT-82/P- 
1434, Jun 1982. 51pp. Presented at Symp. IUTAM-ISSM Modern Developments in Analytical 
Mechanics', Turin, Italy, Jun 1982. 

[17] L. M. Bates, "Geometric quantization of a perturbed Kepler problem," Rept. Math. Phys. 28 
(1989) 289. 

[18] F. Iachello, "Dynamic symmetries and supersymmetries in nuclear physics," Rev. Mod. Phys. 
65 (1993) 569. 



10 



